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ABSTRACT 

Knowledge of the geometry of pulsational modes is a prerequisite for seismic mod¬ 
elling of stars. In the case of slowly pulsating B-type (SPB) pulsators, the simple 
zero-rotation approach so far used for mode identification is usually not valid because 
pulsational frequencies are often of the order of the rotational frequency. Moreover, 
this approach allows us to determine only the spherical harmonic degree, i, while the 
azimuthal order, m, is beyond its reach. On the other hand, because of the density 
of oscillation spectra of SPB stars, knowledge of m is indispensable if one wants to 
assign the radial order, n, to the observed frequency peaks. Including the effects of ro¬ 
tation via the traditional approximation, we perform identification of the mode angular 
numbers {£, m) for 31 SPB stars with available multicolour time series photometry. 
Simultaneously, constraints on the rotational velocity, Vrot) and the inclination angle, 
i, are determined assuming uniform rotation and a constant value of Vmt sini. Depen¬ 
dence of the results on the adopted model is tested using HD 21071 as an example. 
Despite some model uncertainties and limitations of the method, our studies show the 
correct approach to identifying the low-frequency oscillation modes. 
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1 INTRODUCTION 

Slowly pulsating B-type (SPB) stars are main-sequence ob¬ 
jects pulsating in high-order gravity modes. The observa¬ 
tional properties of the SPB stars were broadly discussed, 
e.g., by Balona (20001 and De Cat (2002 20071. These 


modes are driven by the classical K-mechanism operating in 


the metal opacity bump (Dziembowski, Moskalik & Pamy- 
[atnykh 1993 Gautschy fc Saio||1993 1. There are three basic 


problems which complicate the interpretation of the SPB 
pulsations and make seismic modelling of these pulsators 
particularly difficult. First, their oscillation spectra obtained 
from ground-based observations are usually sparse and lack 
equidistant patterns. Therefore, the only way to identify pul¬ 
sation modes, i.e., to determine the spherical harmonic de¬ 
gree, £, and the azimuthal order, m, is the use of the informa¬ 
tion contained in the light and/or line profile variations. The 
exceptions are the stars observed from space, which have os¬ 
cillation spectra with hundreds of frequency peaks and with 


some clear structures (e. g. 

Degroote et al.||2010 

Balona et 

al.||2011| [McNamara, Jackiewicz fc McKeever||20121 [Papics 

et al.|2012 

2014 

1. 


The most popular tools derived from ground-based ob¬ 
servations and used for mode identification are the am- 
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plitude ratios and phase differences in various photomet¬ 
ric passbands (Balona fc Stobie||1979 Stamford & Watson 


|1981[ |Watson|1988| , because these data are much more eas¬ 
ily obtained than the line profiles. The amplitude ratios and 
phase differences in different passbands were called the pho¬ 
tometric non-adiabatic observables by|Cugier, Dziembowski| 


& Pamyatnykh (19941, who improved this method for /3 


Cephei stars by including non-adiabatic effects in their pul¬ 
sational computations. If all effects of rotation are ignored, 
these observables are independent of the azimuthal order, 
m, and inclination angle, i, and modes with different har¬ 
monic degrees, £, occupy disjoint regions on the amplitude 
ratio versus phase difference diagrams. The results for non¬ 
rotating SPB stellar models were presented for the first time 
by Townsend {20021. 

It should be mentioned that the mode degree, £, and 
azimuthal order, m, can be extracted from spectroscopic 
observations as well. For details on the spectroscopic meth¬ 
ods, we refer to, e. g. Campos & Smith ( 1980a|b I; Balona 


Daszyhska (20011; Gies & Kullavanijaya (19881; Kennelly 


fc Walker (19961; Telling fc Schrijvers (19971; Schrijvers fc 


Telling (1999 20021 and Zima (20061. Some of these meth- 


( 1986a|b I; Aerts, de Pauw fc Waelkens (1992l; Cugier fc 


ods take into account the effects of the Coriolis force (e.g.. 
Telling fc Schrijvers|| 19971 ISchrijvers fc Telting||1999] |2002| 


Zima|2006 1 but only to the first order and therefore they do 

not apply to low-frequency g modes. The traditional approx- 
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imation was applied only in the moment method by |De Cat| 
et al. (20051 to a few SPB stars. The most recent overview 


of the spectroscopic methods can be found in |Uytterhoeven| 
(20141). 


The second problem is associated with a very high den¬ 
sity of the theoretical oscillation spectra of SPB stellar mod¬ 
els. For this reason, the assignment of the radial order, n, to 
the individual observed peaks is usually impossible or can 
lead to the misinterpretation of the data. For example, |Deg-| 
roote et al. (2010 2012 I claimed that they have discovered a 


series of dipole zonal modes equidistant in period in the SPB 
star HD 50230. Our findings, however, contradict this result 


(Szewczuk, Daszynska-Daszkiewicz & Dziembowski 2014). 


We showed that this equidistant pattern is rather acciden¬ 
tal and its interpretation in terms of the asymptotic theory 
is dubious. 


The third problem with the interpretation of pulsations 
in SPB stars is rotation, which significantly changes, both, 
the equilibrium model and pulsational properties. In the case 
of these pulsators, rotation plays an especially important 
role because even at slow rotation pulsational frequencies 
can be of the order of the rotational frequency and the effect 
of the Coriolis force cannot be neglected. On the other hand, 
the effect of centrifugal deformation on g modes is small 
(Ballot et al.|20l'2l. The influence of rotation on long-period 


g modes was examined relying on the so-called traditional 


approximation (e. g. 

Lee & Saio 

1997 

Townsend 

2003a|b 

2005 

jDaszynska-Daszkiewicz, Dziembowski & Pamyatnykh| 

2007 

) or on the truncated expansion in Legendre functions 


So far, the mode identification for SPB stars has been 
made mostly in the zero-rotation approximation, e.g., by 
|Townsend| ( |2002l ) and |De Cat et aL] ( |2004p , and the results 
indicated that most of these stars pulsate in dipole modes. 
The only SPB star for which the effects of rotation were in¬ 
cluded in the photometric method of mode identification was 


Eri (^Daszynska-Daszkiewicz, Dziembowski & Pamyatnykh 


2008 Daszyhska-Daszkiewicz et al. 20151. This identifica¬ 


tion was based on linear non-adiabatic calculations with the 
effects of rotation treated in the traditional approximation. 


Therefore, there is a need for this type of studies and 
the main goal of the present work is to make mode iden¬ 
tification with the effects of rotation taken into account 
for all SPB stars with available and appropriate photomet¬ 
ric and/or spectroscopic data. Uniform rotation is assumed. 
We follow the method used in [Daszynska-Daszkiewicz et al.| 
( 20081 which, in addition, allows us to constrain also 

the rotational velocity of our programme stars. For compar¬ 
ison, we make also mode identification in the zero-rotation 
approximat ion. 


The structure of the paper is as follows. In Section 2, we 
present observational data of selected SPB stars gathered 
from the literature, and the luminosities and evolutionary 
masses determined by us. In Sections, the methods of mode 
identification we use, dependence on model parameters and 
the results are discussed. We end with a summary in the last 
section. 


2 PROGRAMME STARS 


We selected 31 SPB stars for which the data needed to put 
a star on the H-R diagram and perform mode identifica¬ 
tion were available. The HD numbers of the stars and their 
stellar parameters are given in Table[^ All the parameters, 
except the masses and luminosities, were taken from the lit¬ 
erature. If there were no determinations of the metallicity 


or colour excess, the solar value Zq — 0.0134 (Asplund et 
ar]|2009 1 and 0 mag were assumed, respectively. In the case 


of HD 28114, the colour excess, E{B — V), in the Johnson 
filters was computed from the colour excess, E{b—y), in the 

I. For the star HD 179588, 


Stromgren filters ( Vogt et ah] 

[Grillo et al.| ( |1922[ ) give the value of the interstellar extinc¬ 
tion in the Johnson V filter. Ay- To obtain E{B — V), we 
assumed the standard value of the ratio of total to selective 
extinction, equal to 3.1. 

The values of luminosity were determined in this paper 
using the bolometric corrections from Flower (19961. Stel¬ 
lar masses were derived from evolutionary tracks and they 
usually correspond to the values of luminosity and effec¬ 
tive temperature from the centre of the error box. There 
are three exceptions: HD 74195, HD 85953 and HD 123515. 
The masses of these stars come from models which lie 
on the left bottom edges of the error boxes, so that 
they have the following values of the effective temper¬ 
ature and luminosity: log refr=4.243, log I//L0=3.323 for 
HD 74195, logreff=4.312, logL/L0=3.833 for HD 85953, 
and logTefr=4.081, logL/L0=2.484 for HD 123515. This is 
because the models from the centres of the error boxes for 
these three stars are located outside the main sequence. 
Since g modes are strongly damped after the main-sequence 
evolutionary stage, it is rather improbable that these stars 
have already left the main sequence. 

The evolutionary computations were performed using 
the Warsaw-New Jersey evolutionary code ( e. g. |Pamyat-| 
nykh et al.|1998[ ), with the OP opacity tables ( Seaton|2005 \ 


and the latest heavy element mixture of [Asplund et al.| 
( |2009| ) (hereafter AGSS09). We assumed the initial hydrogen 
abundance X — 0.7, the metallicities determined from ob¬ 
servations (column 7, Table and the equatorial rotational 
velocities equal to the measured projected value (column 9, 
Table [^. Overshooting from convective core was not taken 
into account. 


We assumed the minimum value of the rotational ve¬ 
locity and the zero value for convective core overshooting 
because their real values are not known. It has to be em¬ 
phasized that these parameters (like metallicity, initial abun¬ 
dance of hydrogen, etc.) can affect the position and extent 
of evolutionary tracks on the H-R diagram as well as the 
parameters estimated from the fit. For instance, if we as¬ 
sume the often-used value of overshooting, aov = 0.2 Hp, 
luminosities and effective temperatures from the centres of 
the error boxes of HD 74195, HD 85953 and HD 123515 are 
in agreement with the main-sequence stages of evolution in 
contrast to the case without overshooting. From these mod¬ 
els and for the central values of log Teft and log L/L 0 , we ob¬ 
tained masses of 6.0 M© for HD 74195, 9.1 M 0 for HD 85953 
and 3 . 8 M 0 for HD 123515 compared to the values 6.1 M 0 , 
8.8 M 0 and 3.7 M 0 from the left bottom edges with ne¬ 
glected overshooting, respectively. The models with over¬ 
shooting give also different maximum values of the rotational 
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Modes identification in rotating SPB stars 3 


Table 1. The astrophysical parameters of the 31 selected SPB stars. Brightness in the second column is given in the Johnson V passband. 
The parameters in the next eight columns have their usual meanings. Remarks are given in the last column. 


HD 

V 

SpT 

logTeff 

logi/L© 

TT 

Z 

E(B-V) 

Vrot sin i 

M 

rmk 


(mag) 


(dex) 

(dex) 

(mas) 

(dex) 

(mag) 

(km s“®) 

(M©) 


1976 

5.581 

B5IV1 

4.199® 

3.295 

3.26® 

0.0165® 

0.091® 

143'® 

6.2 

SBl®® 




±0.005 

±0.176 

±0.63 

+0.0008 

-0.0007 

±0.004 

±9 


VB®° 

3379 

5.871 

B2.5IV1 

4.3011® 

3.397 

3.14® 

0.0134* 

0.070®® 

48®® 

7.1 




B3Vi 

±0.022 

±0.101 

±0.30 



±8 



21071 

6.081 

B5IV1 

4.164® 

2.444 

6.27® 

0.0082® 

0.085® 

50®'® 

3.7 




B7Vi 

±0.007 

±0.076 

±0.40 

+0.0016 

-0.0014 

±0.005 




24587 

4.64‘i 

B5V® 

4 .I 6 II® 

2.626 

8.65® 

0.0134* 

0* 

30®^ 

4.2 

SBl®® 




TO.030 

±0.060 

±0.30 






25558 

5.33® 

B5V'^ 

4.235® 

3.117 

5.09® 

0.0096® 

0.112® 

28®® 

5.5 





±0.008 

±0.071 

±0.32 

+0.0019 

-0.0016 

±0.010 

±2 



26326 

5.43® 

B5IV5 

4.176® 

2.824 

5.00® 

0.0103® 

0.004® 

17®® 

4.5 





±0.006 

±0.070 

±0.28 

+0.0018 

-0.0015 

±0.004 

±1 



28114 

6.061 

B 6 IVI 

4.1461® 

2.672 

5.05® 

0.0134* 

0.149®® 

21®® 

4.3 





±0.031 

±0.099 

±0.49 



±4 



28475 

6.79® 

B5V10 

4.188® 

2.733 

4.31® 

0.0134* 

0.240® 

30®9 

4.6 

SB2®5 




±0.012 

±0.129 

±0.58 


±0.010 




34798 

6.361 

B5IV/V1 

4.197® 

2.766 

3.80® 

0.0115® 

0.038® 

30®^ 

4.6 

VB®° 



B2IV/V1 

±0.006 

±0.119 

±0.47 

+0.0020 

-0.0017 

±0.004 

±9 





B3Vi 











B5Vi 









37151 

7.3811 

B8Vi® 

4.108® 

1.895 

5.20® 

0.0103® 

0.037® 

20®® 

2.8 



7.381® 

B8i® 

±0.008 

±0.109 

±0.62 

+0.0010 

-0.0009 

±0.005 

±2 



45284 

4.1761® 

2.476 

3.02® 

0.0134* 

0* 

52®® 

4.1 





±0.014 

±0.166 

±0.55 



±5 



53921 

5.611 

B9IV1 

4.I37I® 

2.388 

6.61® 

0.0134* 

0* 

17®® 

3.7 

VB®s 



B9III1 

±0.002 

±0.056 

±0.32 



±10 



74195 

3.621 

B3/5V1 

4.217I® 

3.392 

6.61® 

0.0100® 

0.006® 

18®® 

6.1 




BSiyi 

B3III1 

±0.026 

±0.069 

±0.35 

±0.0007 

±0.004 

±2 



74560 

4.841 

Bsiyi 

4.227® 

2.919 

6.73® 

0.0177® 

0.004® 

45®° 

5.4 

SB1®5 



B3Vi 

B4IV1 

B5IV1 

±0.006 

±0.048 

±0.17 

+0.0030 

-0.0026 

±0.004 



VB®° 

85953 

5.93I 

B2III1 

4.3OII® 

3.973 

1.42 

0.0134* 

0* 

30®2 

00 

bo 




B2/3III1 

±0.011 

±0.140 

±0.21 



±1 





B2Vi 









92287 

5.901 

Bsiyi 

4.215I® 

3.178 

3.24® 

0.0134* 

0.070®® 

41®° 

5.8 

SB1®5 



B3III1 

±0.004 

±0.076 

±0.21 



±16 



121190 

5.65®'^ 

B8V®® 

4.0824® 

1.992 

9.24® 

0.0134* 

0.010®® 

118'®° 

3.0 




B9V®® 

±0.010 

0.040 

±0.28 



±3 



123515 

5.981 

B8Vi 

4.079I® 

2.574 

3.98® 

0.0134* 

0* 

24®® 

3.7 

SB2®5 



B 8 IVI 

B9IV1 

±0.002 

±0.090 

±0.39 





VB®° 

138764 

5.161 

B8Vi 

4.I6II® 

2.541 

8.17® 

0.0134* 

0.060®® 

21®® 

4.1 




B 6 IVI 

B7IV1 

±0.030 

±0.081 

±0.58 



±2 



140873 

5.39I 

B 8 IV/VI 

4,14419 

2.583 

7.25® 

0.0134* 

0.116®'® 

80®° 

4.1 

SB2®® 



B 8 IVI 

B 8 IIII 

±0.003 

±0.063 

±0.31 





VB®° 

143309 

9.2814 

B8Ib/Ili5 

4.147®® 

2.557 

1.30®! 

0.0134* 

0.140®® 

10®® 

4.1 




B9Ib/Ili® 

±0.020 

±0.052 







160124 

7.171 

B3IV/V1 

4.171®® 

3.026 

2.16® 

0.0134* 

0.150®®’ 

820 

5.2 




Bsiyi 

±0.020 

±0.239 

±0.58 



±4 





B5Vi 









160762 

3.8OI 

BSiyi 

4.2901® 

3.396 

7.17® 

0.0134* 

0.000®® 

6®4 

6.9 

SBl®® 



B3Vi 

B3.5III1 

±0.011 

±0.059 

±0.13 



±1 


VB®° 

177863 

6.301 

B8Vi 

4.153® 

2.520 

5.38® 

0.0132® 

0.126® 

60®7 

4.0 

SBl®® 



00 

±0.012 

±0.137 

±0.75 

+0.0030 

-0.0024 

±0.029 



VB®° 
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Table 1 — continued 


HD 

V 

SpT 

log Tgff 

logL/L© 

TT 

Z 

E{B - V) 

Trot sin i 

M 

rmk 


(mag) 


(dex) 

(dex) 

(mas) 

(dex) 

(mag) 

(km s“l) 

(M©) 


179588 

6,75^ 

B9IV1 

4.04818 

2.307 

4.12® 

O.OI 34 I 8 

0.02928 

3517 

3.4 

SB 18 ® 



B 8 VI 

±0.009 

±0.180 

±0.84 

+0.0034 

-0.0027 




VB86 

181558 

6.26® 

B5III5 

4.1742 

2.749 

4.178 

0.01352 

0.0772 

2017 

4.6 

VB86 




±0.006 

±0.094 

±0.38 

+0.0027 

-0.0022 

±0.005 




182255 

5.19^ 

B 6 IIP 

4.1482 

2.414 

8.308 

0.01002 

0.0022 

2517 

3.7 





±0.004 

±0.081 

±0.60 

+0.0017 

-0.0015 

±0.001 

±9 



191295 

7.168 

B7III8 

4.11412 

2.435 

2.988 

0.0134* 

0* 

1612 

3.7 





±0.033 

±0.141 

±0.47 



±1 



206540 

6 .O 51 

B5IV1 

4.13012 

2.586 

4.668 

0.0134* 

0.04824 

1512 

4.1 




B7IIP 

±0.032 

±0.078 

±0.37 



±2 



208057 

5.08* 

B3V8 

4.2292 

3.070 

5.168 

0.01102 

0.0202 

11017 

5.5 





±0.005 

±0.058 

±0.24 

±0.0005 

±0.004 

±16 



215573 

5 . 34 I 

B5/7V1 

4.1532 

2.591 

6.628 

0.01102 

0.0172 

812 

4.0 




Beiv^ 

±0.003 

±0.057 

±0.19 

+0.0008 

-0.0007 

±0.003 

±1 




(i)|R ee d|(|2005^; (2)|Niemczura|l|2003[l; ( 3) [van Leeuweii|(|2007^ ; (4)|Huang, Gies &: McSwain|)|2010[l; (4)|CIausen et al.|([T997^; (5)|Houk | 
I&: Smith-Moor^^l988|l; (6) |DucaVi| ( |2002^; [Cucchiaro et aT!| *| |1980|l; (8) |Hohle, Neuhauser fc Schutz| l|2010[); (9) |Fitzpatrick fc Massa 
( |2007[l; (10)|Cowley et alT]||1969[ l; (ll)|Warren fc Hesser|||1977[ l; ( 12)|Lefever et al.|||2010[ |; (13) |Koen &: Eyer|||2002|; ( 14) |Kilambi|||1975[l; 


(15)|Houk Sz Cowley|(|1975|l; (16) |Gudennavar et aL|||2012|); (17)|Abt, Levato Sz Grosso|(|2002||; (18)|Wu et al.|( 

2011|l; (19) {Briquet et al.| 

<|2007||; (20)|Hubrig et al.|(|2006|l; (21) |Kharchenko et al.|' 

2009f; (22)|Vogt et al.|(|1998|l; (23)|GriUo et al.|(|1922 

l; (24) 

Tobin 

(|1985|l; (25) 

TurnerHl973|l; (26) |Wegner| (|2003|l; (27) |Vleeming| (|1974 

i; (28)|Nieva & Przybilla| (|2012||; (29) |Strom, Wolff & Dror 

■poos 

; (30)|Aerts| 

et al.|(|1999|1; (31)|Hubrig et al.| 1|2009||; (32)|Zorec Sz Royer|(|2012||; (33)|Grosso Sz Levato|||2011|1; (34)|Bailey Sz Landstreet 

o 

CO 

5 (35) 

Pourbaix et al.|||2004|l; (36)|Mason et al.|(|2001||; (37)|Slawson Sz Landstreet] ||1992|); (38)|Whittet Sz van Breda]||1980|l; (39)|Houk|i 

19781); 

(40) Aerts & Kolenberg i 2005 i; (*) Zq = 0.0134 assumed. 


frequencies, r'jot* = 0.7623 d ^ for HD 74195, = 0.6147 

for HD 85953 and = 0.9576 d“^ for HD 123515 
compared to = 1.0388 d“^, = 0.8281 d“^ and 

I'rot* = 1.1284 d~^ obtained from models from left-bottom 
edges with neglected overshooting, respectively. Therefore, 
based on the differences between the cases with and without 
overshooting, we can estimate the error in mass in Table[^ 
to be at least a few tenths of a solar mass, and the error in 
the critical rotation frequency in Tablej^to be at least 0.25 
d“^. 

In Fig. 1, the programme stars are plotted in the H- 
R diagram. Also shown is the SPB instability strip for the 
dipole and quadrupole modes. To determine the extent and 
position of the instability strip, we used models with the 
parameters as described earlier, but with the metallicity 
Z = 0.015 and with the rotational velocity Vrot=0 km s“^. 
The numbers are assigned to the individual stars. It should 
be noted that due to various metallicities and rotational 
velocities, the models used to estimate masses listed in Ta- 
ble[3 do not correspond strictly to the evolutionary tracks 
presented in Fig. 1. 


3 MODE IDENTIFICATION 

In order to make mode identification, we will use observa¬ 
tions in the UBV Geneva filters. We include only these 
three filters because the response functions in the Geneva 
photometric system widely overlap and all seven filters do 
not give independent information. Therefore, we chose the 
UBV filters which have non-overlapping response functions. 
In Table|^ we list all frequencies detected in the programme 


stars, and the corresponding values of the amplitudes in the 
V filter. In the two second last columns, we give also the 
minimum value of the rotational frequency, resulting from 
Kotsini, and the critical rotational frequency correspond¬ 
ing to the break-up velocity values of a 

mass M and a radius R were adopted for the central models 
of Table 1. 

In this section, we present the results of mode identih- 
cation for the frequency peaks observed in the selected SPB 
stars. First, we consider the case when all effects of rotation 
are neglected. Then, we give results in the framework of the 
traditional approximation in which the effect of the Coriolis 
force on long-period g modes is included. 


3.1 Methods 


3.1.1 The zero-rotation approximation 


The most common way to identify pulsation modes from 
photometry is to compare the empirical amplitude ratios and 
phase differences in different photometric passbands with 
their theoretical counterparts. Usually, we do this by com¬ 
puting a discriminant (e. g. |Daszyhska-Daszkiewicz fc Wal-| 
cza ^[2009t 


2 


X 


1 

2(iV- 1) 


2(iV-l) 

E ' 


{X° 


i=l 



( 1 ) 


where Af and Xj are the observational and theoretical pho¬ 
tometric observables, respectively, i.e., the amplitude ratios, 
AxjAy, or phase differences, cjix — (py: in V photometric 
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'4.19 4.18 4.17 4.16 4.15 4.14 4.13 


log Teff 


Figure 1. The programme stars (cf.TableQ in the H-R diagram. A zoom-in on the central, densely populated area around logTgff = 4.16 
and logL/LQ = 2.6 is shown in the bottom panel for clarity. The SPB instability strip for the dipole and quadrupole modes is shown as 
a grey area. Stability calculations were performed with the initial abundance of hydrogen X = 0.7, metallicity Z = 0.015, OP opacities, 
AGSS09 element mixture, for non-rotating models, and without taking into account overshooting from the convective core. There are 
shown the evolutionary tracks for the masses 3Mq (solid line), 5Mq (dashed line) and 7Mq (dotted line). The labels (1—31) indicate 
the following stars: 1 - HD 1976; 2 - HD 3379; 3 - HD 21071; 4 - HD 24587; 5 - HD 25558; 6 - HD 26326; 7 - HD 28114; 8 - HD 28475; 9 - 
HD 34798; 10 - HD 37151; 11 - HD 45284; 12 - HD 53921; 13 - HD 74195; 14 - HD 74560; 15 - HD 85953; 16 - HD 92287; 17 - HD 121190; 
18 - HD 123515; 19 - HD 138764; 20 - HD 140873; 21 - HD 143309; 22 - HD 160124; 23 - HD 160762; 24 - HD 177863; 25 - HD 179588; 26 
- HD 181558; 27 - HD 182255; 28 - HD 191295; 29 - HD 206540; 30 - HD 208057; 31 - HD 215573. 
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6 W. Szewczuk and J. Daszynska-Daszkiewicz 


Table 2. The frequencies and amplitudes in the V Geneva filter of the selected SPB stars. Symbols p and s denote that the frequency is 
detected in the photometric and spectroscopic observations, respectively. The two second last columns contain the minimum and critical 
values of the rotational frequency. 


- (d-i) 


HD 


Av 

(mag) 


(d-i) 

(d-i) 

Ref 

1976 

0.39946(4)P 

0.93895(4)P 

1.20346(6)P 


0.4770 

0.8138 

P- (1) 


6.9(6) 

5.9(6) 

4.0(6) 





3379 

1.82023(19)P 

1.59418(22)P 



0.2273 

1.4643 

P- (1) 


3.8(6) 

3.6(6) 






21071 

1.18843(1)P= 

1.14934(3)P 

1.41968(7)P 

0.95706(9)P 

0.3768 

2.1256 

P- (1) 


18.5(6) 

7.7(6) 

3.8(6) 

3.0(6) 



s- (2) 

24587 

1.1569(6)P= 




0.1814 

1.6385 

ps - (3) 


6.6(8) 







25558 

0.65265(2)P= 

1.93235(8)P 

1.17913(10)P 


0.1351 

1.3287 

P- (1) 


14.7(6) 

3.8(7) 

3.6(6) 




s- (2) 

26326 

0.5338(8)P= 

0.1723(8)'* 

0.7629(10)P 


0.0875 

1.3296 

ps - (3) 


8.3(1.0) 

- 

7.0(1.0) 





28114 

0.48790(4)P 

0.48666(8)P 



0.1123 

1.3671 

P- (1) 


14.4(9) 

7.2(1.0) 






28475 

0.68369(7)P 

0.40893(13)P 



0.1810 

1.6978 

P- (1) 


8.4(1.0) 

5.3(1.0) 






34798 

0.78350P 

0.85231P 

0.97432P 

0.67795P 

0.1818 

1.7114 

P - (2) 


12.6(7) 

2.4(7) 

5.9(7) 

7.0(7) 




37151 

1.243166P 

1.180186P 

1.042663P 

1.105252P 

0.2195 

3.2770 

P- (4) 


9.5(4) 

1.452063P 

9.8(4) 

8.7(4) 

7.7(4) 





1.9(4) 







45284 

1.23852P 

1.12753P 

1.50605P 


0.4006 

2.3120 

P - (2) 


7.4(7) 

10.5(8) 

4.7(7) 





53921 

0.6054(6)P= 




0.1209 

1.9618 

ps - (3) 


5.3(7) 







74195 

0.35475(9)P= 

0.35033(9)P= 

0.34630(9)P'* 

0.39864(9)P** 

0.0710 

1.0388 

ps - (3) 


13.2(4) 

10.6(4) 

3.5(4) 

3.0(4) 




74560 

0.644714(2)P= 

0.395771(10)P 

0.447646(21)P 

0.635982(33)P 

0.2632 

1.7639 

ps - (*) 


15.0(4) 

4.2(4) 

3.3(4) 

2.2(4) 





0.822845(43)P 

1.610567(54)P 







3.1(4) 

2.6(4) 






85953 

0.2663(6)P= 

0.2189(6)P 

0.2353(7)** 


0.0903 

0.8281 

ps - (3) 


8.5(7) 

3.8(7) 

- 





92287 

0.21480(7)P“ 




0.1681 

1.0669 

ps - (3) 


10.2(6) 







121190 

2.6831(4)P= 

2.6199(4)P 

2.4713(7)P 


1.0273 

2.3831 

ps - (4) 


4.8(7) 

3.9(6) 

2.6(8) 





123515 

0.68528(10)P= 

0.65929(10)P 

0.72585(10)P** 

0.55198(10)P 

0.1181 

1.1284 

ps - (3) 


16.4(5) 

11.3(5) 

9.3(5) 

5.2(5) 




138764 

0.7944(9)P= 

0.6372(9)** 



0.1395 

1.8539 

ps - (3) 


17.7(1.1) 

- 






140873 

1.1515(8)P= 




0.4690 

1.5388 

ps - (3) 


15.4(1.6) 







143309 

0.59830(2)P 

0.60213(2)P 

0.71414(2)P 

1.20066(2)P 

0.0611 

1.6277 

P - (2) 


23.0(1.5) 

0.59133(2)P 

17.6(1.5) 

10.7(1.5) 

6.7(1.4) 





7.6{1.4) 







160124 

0.52014P 

0.52096P 

0.52055P 

0.52137P 

0.0319 

0.9692 

P - (2) 


11.6(1.6) 

10.1(1.4) 

12.7(1.4) 

11.0(1.5) 





0.52197P 

0.70259P 

1.0422P 

0.31366P 





12.3(1.5) 

7.6(1.3) 

5.8(1.3) 

5.7(1.3) 




160762 

0.28671P 




0.0360 

1.3462 

P - (2) 


7.8(2.9) 







177863 

0.84059(10)P= 

0.10108(10)P 



0.3943 

1.8074 

ps - (3) 


17.4(4) 

2.8(4) 






179588 

0.85654(4)P 

2.04263(5)P 

2.19989(7)P 

1.83359(9)P 

0.1815 

1.1582 

P- (1) 


12.9(1.3) 

10.5(1.3) 

8.0(1.3) 

5.9(1.4) 




181558 

0.80780(10)P= 




0.1114 

1.4997 

ps - (3) 


27.5(7) 
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Table 2 — continued 




u (d- 






HD 


Av (mag) 


(d-i) 

(d-fi 

Ref 

182255 

0.97185(3)P= 

0.79225(4)P= 

0.62526(6)P 

1.12780(10)P 

0.1815 

2.0047 

P - (1) 


18.7(9) 

1.02884(11)P 

12.7(9) 

8.3(8) 

5.8(8) 



s- (2) 


4.3(1.1) 







191295 

0.71908(6)P 

0.49011(12)P 

0.29302(17)P 


0.0972 

1.5462 

P - (1) 


15.7(9) 

8.6(1.0) 

5.8(9) 





206540 

0.72002(4)P 

0.62125(5)P 

0.38271(8)P 


0.0822 

1.3811 

P - (1) 


12.2(1.1) 

8.6(1.1) 

5.6(1.0) 





208057 

0.89050(7)P= 

0.80213(10)P= 

2.47585(8)P 


0.5440 

1.3741 

P - (1) 


8.4(1.4) 

6.0(1.3) 

7.4(1.5) 




s- (2) 

215573 

0.5654(6)P'’ 




0.0484 

1.5973 

ps - (3) 


14.9(2.0) 








(i)|^ Cat et ar] | |2007| l; (2) |De Cat| | |20d2l l; (3) |De Cat &c Aerts| ) [^02| l; (4) [North fc Paltani| | [T994[ | ; (4) [Aerts &: Kolenbergj l |2005| l; 
(*) this paper. 


passbands. N denotes the number of passbands and ai are 
the observational errors of the quantities X°. 

If all effects of rotation on stellar pulsation are ne¬ 
glected, the photometric observables are independent of the 
inclination angle, i, intrinsic mode amplitude, e, and the 
azimuthal order, m. Thus, the method allows us to deter¬ 
mine only the mode degree. This is the main disadvantage 
of disregarding the rotational effects. 

In the mode identification process, we accept the degree 
£ if the corresponding value of is below some specified 
level. We choose this level in such a way that the probability 
of rejecting a correct solution is less than 20 per cent. In 
the case of three passbands, it is about ~ 1-5- We note, 
however, that it is quite a rough estimate, because we do 
not deal here with a strict statistical approach. Moreover, 
due to an increase of the averaging over the stellar disc with 
the increasing mode degree, we considered modes with I up 
to 6. 


The formulae for the complex amplitude of the light 
variation in the x passband in the framework of the linear 
theory in the zero-rotation approximation were given, e. g. 


by Daszyhska-Daszkiewicz et al. (20021 and are not repeated 
here. 

To compute the theoretical values of the photometric 
amplitudes and phases, input from models of stellar at¬ 
mospheres and linear non-adiabatic pulsation calculations 
are needed. The former include the flux derivatives over 
effective temperature and gravity and the limb-darkening 
coefficients and their derivatives. While the latter, the so- 
called non-adiabatic parameter / dehned as the ratio of 
the amplitude of the radiative flux perturbation to the ra¬ 
dial displacement at the photosphere level. For details see, 
e. g. jPaszynska-Daszkiewicz, Dziembowski fc Pamyatnykh] 


(20031 and Daszyhska-Daszkiewicz & Szewczuk ( 2011| . In 
what follows, we rely on the Kurucz (20041 stellar atmo¬ 
sphere models. 

A significant reduction in the number of modes to be 
considered comes from information about their instabil¬ 
ity. First, we took into account only theoretically unstable 
modes, i. e., those with the instability parameter Tj > 0. 
However, because of uncertainties in the model and opac¬ 


ity data, this condition was slightly relaxed, and ultimately 
we accepted the modes with rj > —0.05. 


3.1.2 Including rotation 


Taking into account rotation in the mode identihcation de¬ 
mands a more sophisticated treatment than that described 
in the preceding subsection. We shall use the method pro¬ 
posed by Daszyhska-Daszkiewicz et al. (2008 20151, which 
includes the effects of rotation via the traditional approxi¬ 
mation. In this method, the observed complex amplitudes 
are compared with the corresponding theoretical values us¬ 
ing the statistic: 


1 ^ 

where Ai and Ti are the observed and theoretical complex 
amplitudes, respectively. The values of Wi = l/uf are the 
statistical weights determined by the observational errors, 
(Ji. N denotes the number of passbands and M is the number 
of unknowns to be determined. The value of the intrinsic 
mode amplitude, e, can be obtained from observations by 
minimizing 


^ E!Ii mArT* 

I T* 12 ’ 

The Hnal expression for the mode discriminant is 
ten as 


(3) 


writ- 


2 


X 


1 

2N- M 




E N I x* 12 I ’ 


where the asterisk denotes complex conjugate. 

Now, the value of x^ depends on both mode numbers {£, 
m), the inclination angle i and rotational velocity Kot- To 
reduce the number of unknowns, the projected rotational 
velocity is kept constant i. e., Kotsini = const. Thus, in 
principle, the method allows deriving simultaneously a value 
of the rotational velocity, Kot- 
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5 

4 

3 

2 

1 

0 


log Te„=4.164 log \JLr,=2AA 
log Tg„=4.171 log L/Lg=2.37 {Model#1) 
log Te„=4.157 log L/Lg=2.37 {Model#2) 
log Tg„=4.171 log L/Lq= 2.52 {Model#3) 
log Tg„=4.157 log L/Lq= 2.52 {Model#4) 


v-|=1.18843 [d‘^] 


0 

-1 




1 2 3 , 4 5 6 


Figure 2. The discriminant ^ function of the mode degree 

£, for the dominant frequency of HD 21071, calculated in the zero- 
rotation approximation for models from the centre (the solid line) 
and four edges (the other lines) of the error box. The black hor¬ 
izontal line indicates our confidence level (x^=1.5). In the lower 
panel, there is shown the corresponding value of the instability 
parameter, rj. 


The theoretical amplitude Jq is calculated with the ef¬ 
fects of rotation taken into account. The expression for 


can be found in Daszynska-Daszkiewicz et al. (2007 20151 


and need not be repeated here. Similarly to the zero-rotation 
case, two inputs are needed. We used the same atmosphere 
models as before, but pulsational calculations which give us 
the /-parameter are performed in the framework of the tra¬ 
ditional approximation in which Coriolis force is included. 

One of the assumptions of the traditional approxima¬ 
tion is that centrifugal distortion of a star is negligible. 
Therefore, we considered the angular frequency of rota- 


condition corresponds approximately to the limit of the half 
of the critical equatorial velocity, (assuming that at 

Q = Dcrit we have Rpoie = 2/37?eq, where 7?poiB and R^q are 
the pole and equatorial radius, respectively). 

As previously, only modes with ^ ^ 6 are assumed to 
be visible in the ground-based photometry. However, this 
time we are dealing with the rotational splitting 2£+l, so 
we have to consider all modes with \m\ ^ 1. In addition, 
we include mixed-Rossby modes, which may be excited and 
visible in the light variations once rotation is fast enough. 
They are retrograde with m = —I and are sometimes called r 
modes ( |Lee|2006j|Daszyhska-Daszkiewicz et al.|2007| ). Again 
only modes with r) > —0.05 were considered. The maximum 
accepted value of the is specified in the same way as in 
the case without rotation. Moreover, |Daszyhska-Daszkiewic^ 


et al. (2008 20151 showed that intrinsic mode amplitude in 


the SPB stars should not be greater than 0.01. Therefore, 
we rejected all identifications with |£| > 0 . 01 . 

If a given mode is also seen in the radial velocity varia¬ 
tions, one can use the information about the mode geometry 
contained in the amplitude of this variations. Only modes 
which reproduce the amplitude of the radial velocity varia¬ 
tions, within the observational error are considered as 

responsible for a given frequency peak. The value of 
was calculated with the empirical value of e obtained from 
photometry (see equation 3). In the case of a star in which 


some frequencies are visible in, both, photometry and spec¬ 
troscopy and some only in the photometry, we adopted a 
fairly reasonable assumption that the radial velocity ampli¬ 
tudes of modes invisible in spectroscopy should be smaller 
than the smallest observed amplitude. Otherwise, such fre¬ 
quencies would be detected in spectroscopy. 


3.2 Working example: HD 21071 

Our working example is HD 21071 (V576 Per), the SPB star 
with four frequency peaks observed in photometry (see Ta¬ 
ble]^. The dominant frequency was also detected in the 
spectroscopic data. As for the other stars, identification of 
pulsation modes in the case of HD 21071 was made using 
the evolutionary models from which the masses given in the 
10th column of Table [T]were estimated. 

We started our analysis with the zero-rotation approx¬ 
imation. In Fig.[^ we show the value of the discriminant 
for the dominant frequency observed in HD 21071 as a func¬ 
tion of the mode degree, £. The models corresponding to the 
centre and four edges of the error box were considered. As 
one can see, a deep minimum at f = 1 is present regard¬ 
less of the adopted model parameters. Thus, the frequency 
zzi is undoubtedly that of a dipole mode. The value of the 
instability parameter, 77 , for the £ = 1 mode is close to 0 
and therefore acceptable. The results of mode identification 
in the zero-rotation approximation for other frequencies are 
given in the second column of Table The results for the 
other frequencies are usually model independent with the 
exception of U 2 and 1/3 for which we obtained £ = 1 or 2 if pa¬ 
rameters of Model#l or Model 7(^:2 were assumed. The effect 
of core overshooting on identification of £ was tested using 
Model^^tS^ i.e., the central model with Oov = 0.2 Hp. This re¬ 
sult is not shown in Fig. [^because for fhi® model almost 

coincides with for the central model without overshoot¬ 
ing and they are hardly distinguishable. Results for other 
frequencies are also compatible with those from the central 
model without core overshooting. 

Mode identification with the effects of rotation taken 
into account consists of several steps. Calculations were per¬ 
formed in the range of rotational velocities from 50 to 250 
km s“^ corresponding to the observed value of Kot sini and 
about a half of the critical rotational velocity for HD 21071, 
respectively. First, we rejected all modes with rj < —0.05. 
In Fig.[^ we show an example of the instability parameter 
for the dipole, quadrupole and r modes with m = — 1 , -2 
for the three values of the rotational velocity, Kot = 50, 
150, 250 km s“^, as well as for the zero-rotation case. A 
model of HD 21071 with the central values of the error box, 
logTefi = 4.164, logL/L 0 = 2.444, M = 3.69Mq, was as¬ 
sumed. In Fig. 3, we mark also frequencies and the Geneva 
V amplitudes of the four peaks detected in HD 21071. Rely¬ 
ing only on the instability information, one can see that the 
observed frequencies can be associated with some modes at 
certain values of the rotational velocities. 

The frequency vq can be associated with the modes (1, 
0 ) or ( 2 ,- 1 ) in the whole range of the rotational velocities. 
The modes (1, -|-1) or (2, 0) are possible only for the lowest 
values of Kot, and (r, — 2 ) for a narrow range of K-ot around 
120 km s“^. The peak vq could be also the (2, —2) mode, 
but this case is quite tricky. In the lower limit of the rotation 
rate, the (2, —2) mode is unstable at the frequency v\. When 
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V4V2V-, V3 


V4V2V1 V3 



0 [km s'^] 
50 [km s'^] 
150 [km s'^] 
250 [km s"''] 


0 0.5 1 1.5 2 2.5 3 

V [d'l 


Figure 3. The instability parameter for the dipole, quadrupole and r modes with m = —1, -2 in the central model of HD 21071 
(logTeff = 4.164; logL/LQ = 2.444; M = 3.69Mq). Four values of the rotational velocity were considered: 0 km s~^ (open diamonds), 
50 km s“^ (filled squares), 150 km s~^ (filled circles) and 250 km s“^ (filled triangles). The observed Geneva ^-filter amplitudes are 
plotted with the vertical solid lines. 
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inclination [°] 



50 100 150 200 250 

Vrot [km s'^] 


Figure 4. Top panel: the value of a function of the rotational velocity for the dominant frequency 18843 of HD 21071, 

calculated using the traditional approximation. The stellar parameters of the model from the centre of the error box were assumed. 
Horizontal dotted line indicates the 80 per cent confidence level. 

Middle panel: the value of the intrinsic mode amplitude, e, as a function of the rotational velocity for the dominant frequency of HD 21071. 
Only modes presented in the top panel are shown. The horizontal dotted line marks the adopted maximum of |£|. 

Bottom panel: the theoretical value of the radial velocity amplitude as a function of the rotational velocity for the dominant frequency 
observed in HD 21071. Only modes presented in the top panel are shown. The observed range of the radial velocity amplitude for vi is 
marked with the horizontal lines. 
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Table 3. The results of mode identification of the frequency peaks observed in HD 21071 from the four approaches described in the text. 
The adopted model parameters are: logTeff = 4.164, logL/LQ = 2.444, M = 3.69M0. The values of the rotational velocity and the 
rotational frequency are in km s“^ and d“^, respectively. Column a — solutions which satisfy the condition; column b — solutions 
which satisfy additionally the e condition; column c — solutions which satisfy the conditions and Ayr', column d — solutions which 

satisfy the conditions x^i Ayr, and the condition on the common range of Vrot- 


u (d ^- 01=0 a b c 


1.18843 1 (1,0) [50;250] (1,0) [55;250] (1,0) [211;250] (1,0) [211;212] {1.5900;1.5976} 

(1,1) [50;90] (1,1) [50;90] 

(2,0) [62;63] 

(3,-1) [55;58] 

(4,-2) [55;55] 


1.14934 1 

(1,0) [50;235] 

(1,0) [51;235] 

(1,0) [51;235] 

(1,0) [211;212] {1.5900;1.5976} 


(1,1) [50;85] 

(1,1) [50;85] 

(1,1) [50;85] 



(2,0) [50;60] 

(2,0) [50;59] 

(2,0) [50;59] 



(2,0) [62;75] 
(3,-1) [54;66] 
(4,-2) [54;55] 

(2,0) [64;75] 




1.41968 1 

(1,0) [135;255] 

(1,0) [135;255] 

(1,0) [135;255] (1,0) [211;212] {1.5900;1.5976} 


(1,1) [60;125] 

(1,1) [60;125] 

(1,1) [60;125] 


(2,0) [50;60] 

(2,0) [50;59] 

(2,0) [50;59] 


(2,0) [62;125] 

(2,0) [64;125] 

(2,0) [69;125] 


(2,-1) [50;50] 

(2,1) [50;60] 

(2,1) [50;60] 


(2,1) [50;60] 

(3,0) [62;65] 

(3,0) [62;65] 


(3,0) [50;50] 

(4,-1) [58;66] 

(4,-1) [58;66] 


(3,0) [62;65] 
(3,-1) [52;63] 
(3,-2) [50;50] 
(4,-1) [50;50] 
(4,-1) [58;81] 
(4,-1) [128;175] 
(4,-2) [54;55] 
(4,-3) [50;50] 
(5,-2) [50;50] 
(5,-2) [62;63] 
(5,-3) [53;53] 
(5,-4) [50;50] 
(6,-5) [50;50] 

(4,-1) [128;175] 

(4,-1) [135;175] 


0.95706 1, 2 

(1,0) [50;155] 

(1,0) [50;155] 

(1,0) [50;155] (2,-1) [211;212] {1.5900;1.5976} 


(1,-1) [50;105] 

(1,-1)[50;105] 

(1,-1) [50;105] 


(1,1) [50;60] 

(1,1) [50;60] 

(1,1) [50;60] 


(2,-1) [50;212] 

(2,-1) [50;212] 

(2,-1) [50;212] 


(2,-2) [50;57] 

(2,-2) [50;57] 

(2,-2) [50;56] 


(2,-2) [180;196] 
(3,-2) [50;50] 
(4,-2) [55;56] 
(4,-3) [50;50] 
(5,-4) [50;50] 
(6,-5) [50;50i 

(2,-2) [180;189] 

(2,-2) [180;189] 


ie,m) [ 1 /” 


d 


the rotation increases, the instability domain shifts towards 
lower frequencies and the mode (2, —2) becomes stable. Fi¬ 
nally, for sufficiently high values of Kot, the mode (2, —2) is 
reflected and its instability domain again encompasses the 
frequency v,. The frequency V 2 is close to v, and therefore its 
instability is approximately the same. The third frequency, 
1 ^ 3 , is the highest. It can be associated with the modes (1, 0), 
(2, —2) for the highest rotation rates, with (1, -1-1), (2, —1), 
(2, 0), (r, —2) for the slow and moderate rotation rates, and 
with (2, -1-1) for the slowest rotation. The lowest observed 
frequency is 124 . In the regime of the low rotation velocity, it 


can be the (1, —1) or (1, -|-1) mode. In the case of the low 
to moderate rotation rate it can be the dipole axisymmetric 
or the (r, —2) mode. 124 can also be the (2, —2) mode for the 
lowest and high rotation rates and the (2, —1) mode in the 
whole range of the Kot- 

The (2, -1-2) mode cannot account for any of the ob¬ 
served frequencies because its instability occurs at higher 
frequencies than the observed ones. None of the peaks may 
be associated with the (r, —1) mode. It is worth to em¬ 
phasize that r modes are unstable only in narrow ranges of 
frequencies. 
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In the next step, we computed the values of using 
equation (4). As an example, we present results for the dom¬ 
inant frequency (see Fig.[^ top panel). Five unstable modes 
met our x^ < 1-5 condition. The angular numbers (^, m) 
of these modes and the ranges of Kot are listed in the third 
column of Table]^ We called this approach the a case. 

In the b case, we rejected modes with the empirical 
values of |£| above the adopted maximum of 0.01. As one can 
see from Fig.[^(the top and middle panels), in the ranges of 
the rotational velocity in which the value of x^ is acceptable, 
the modes (2,0), (3, —1) and (4, —2) have too large values 
of |e|. The modes with the appropriate values of x^ and |£| 
are (1, 0), (1, -1-1); they are listed in the fourth column of 
TableEl 

Another constraint used by us is the value of the ampli¬ 
tude of the radial velocity changes, and this is the c 

case. In the case of HD 21071, only the dominant frequency 
was detected in the spectroscopic data. In Fig.|^(the bottom 
panel) are shown the theoretical values of calculated 

with the values of e obtained from equation (3) for modes 
presented in the top panel. As one can see, the (1, -|-1) mode 
should be rejected due to the excessive value of the Av^^^, 
while the (1, 0) mode has the theoretical value of Ay^^^ 
within the observational error for Vrot above 211 km s“^. 

In the case of the other frequencies, we required that their 
theoretical values of Ay^^^ were smaller than the observed 
Av„d for the dominant frequency. These identifications are 
given in the fifth column of Table[^ 

Finally, we selected only modes with a common range of 
Vrot for all frequencies - the d case. This gives us a very nar¬ 
row range of the rotational velocities, V)ot £ [211; 212] km 
s“^ and unambiguous mode identification for all observed 
frequencies: vi, V 2 , fs are dipole axisymmetric modes and 
Vi is a quadrupole mode with m = — 1 (the last column in 
Table [^. 

To check the dependence of our mode identihcation on 
the adopted model parameters, the above described proce¬ 
dure was repeated for the models corresponding to the four 
edges of the error box and for the model from the centre of 
the error box calculated with overshooting from the convec¬ 
tive core taken into account, ctov = 0.2 Hp. The results are 
presented in Table|^ Depending on the model parameters, 
we obtained the rotation velocity from about 150 to 250 km 
s-i. 

For Model^l and Model#2, we obtained additional so¬ 
lutions for Vi and the smaller values of Kot. In the case of 
Model^3, we were unable to find mode identihcation for Vi 
in the common range of Vrot with other frequencies. There¬ 
fore, in Table|^ are listed in italics the three modes with the 
values of Kot closest to the range of Kot obtained for vi, 

V 2 and v^. The most complicated situation is for the coolest 
and most luminous Model#4. In this model, for vi we were 
able to hnd the modes which ht the amplitude of the radial 
velocity changes only to within 2cr (in Table|^ indicated as 
2 (t) and it gave us two solutions. Another problem is the 
fact that in the case of the solution for vt, which is con¬ 
sistent with other models, namely the (1, 0) mode, there is 
no common range of Vrot with the results for Vi (the three 
modes with the closest rotational velocity ranges are given 
in Tablej^in italic). In the case of the model with core over¬ 
shooting (Model^5), as for Model#3, we were unable to find 


mode identification for Vi in the common range of Kot with 
other frequencies. 

It is worth to mention that in more evolved stars over¬ 
shooting from the convective core may affect a model more 
signihcantly than in the case of HD 21071. Moreover, there 
is a well-known problem in the literature that discrepancies 
between different determinations of the fundamental param¬ 
eters exist (e. g. |Balona et al.|2011[ ). These uncertainties can 
affect mode identification. From Table|^ one can see that 
the results of our mode identification, especially the values 
of the rotational velocity, depend on the adopted model. Un¬ 
certainties in our results are discussed in more details in the 
last section. 


3.3 Other results 


The results of mode identification for the remaining pro¬ 
gramme stars are presented in Table[^ It is clearly seen 
that the best constraints on {i, m) and the rotational ve¬ 
locity were obtained for the stars with the measured radial 
velocity variations. 

Most of the stars we studied in this paper pulsate in the 
normal modes. There are two exceptions: V 2 in HD 28475 and 
1^3 in HD 206540 which may be either an r mode or a normal 
g mode. 

In the case of vt in HD 53921 (identihcation with and 
without rotation), vi in HD 37151 (identihcation without ro¬ 
tation) and V 2 in HD 179588 (identihcation with and without 
rotation), we are not able to hnd solutions with x^ ^1.5 and 
therefore we accepted higher values of x^ (indicated in Table 
[^with an exclamation mark for Xmax=3.3 which corresponds 
to the probability that real solution has higher x^ equal to 
1 per cent or by a double exclamation mark for Xmax=4.6 
which corresponds to the probability that real solution has 
higher x^ equal to 0.1 per cent). We would like to remind 
that these quantities should be treated more like arbitrarily 
chosen thresholds than a strict statistical approach. 

In the case of spectroscopic binaries, it was checked 
whether the pulsation frequency can be a multiple of the 
orbital frequency. Only for the frequency vi in the star 
HD 177863, the ratio of the pulsation frequency to the or¬ 
bital frequency, equal to 10.016±0.001, is close, but still not 
equal to, within observational error, the integer. A possi¬ 
bility that this mode is excited by tidal forces was studied 
by [Willems fc Aerts] ( |2002) ) but they did not conclude which 
mechanism, classical k or tidal forces, is responsible for driv¬ 
ing this mode. Another possibility is that some observed fre¬ 
quencies are the result of rotational modulation due chem¬ 
ical and/or temperature inhomogeneities. The frequencies 
which fall in the rotational frequency range determined by 
us are: vi and V 2 in HD 74560, vi and Vi in HD 123515, V 2 in 
HD 208057, vi and V 2 in HD 28114, vi and V 2 in HD 28475, 
vi and Vi in HD 34798, all frequencies of HD 37151, vi and V 2 
in HD 45284, all frequencies of HD 143309, vi in HD 179588, 
all frequencies of HD 191295 and HD 206540. We presume 
that most of these frequencies are of pulsational origin, but 
it must be remembered that pulsational and rotational vari¬ 
ability can coexist in one star (e. g. Papics et al.|2012 Thoul 
et al.|20l3 1 and even spots alone could cause multi-periodic 


signal if the stellar surface rotates differentially (e. g. Deg- 
roote et al.||2011l. However, to decide with certainty, a de- 
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Table 4. The results of the mode identification of the frequency peaks observed in HD 21071 based on the evolutionary models from the 
four edges of the observational error box and from the centre of the observational error box but with overshooting from convective core 
taken into account. 


t{Vrot=0) {e, m) km s"! 


^ (d-1) 


Model#l 


Model#2 


Model#3 


Model#4 


Model#5 

1.18843 

1 

(1,0) [148;170] 

1 

(1,0) [174;200] 

1 

(1,0) [249;255] 

1 

(1,0) [229,245]2‘^ 
(1,1) [107;110]2'" 

1 

(1,0) [212;235] 

1.14934 

1 

(1,0) [148;170] 

1 

(1,0) [174;200] 

1 

(1,0) [249;255] 

1 

(1,0) [229,245] 

1 

(1,0) [212;235] 


2 


2 





(1,1) [107;110] 



1.41968 

1 

(1,0) [148;170] 

1 

(1,0) [174;200] 

1 

(1,0) [249;255] 

1 

(1,0) [229,245] 

1 

(1,0) [212;235] 


2 







(1,1) [107;110] 
(2,0) [108;110] 



0.95706 

1 

(1,-1) [148;170] 

1 

(1,-1) [174;200] 

1 

(1,0) [50; 180] 

1 

(1,0) [50; 215] 

1 

(1,0) [50;155] 


2 

(2,-1) [148;170] 

2 

(2,-1) [174;200] 

2 

(2,-1) [50;189j 

2 

(2,-1) [50; 177] 

2 

(2,-1) [50;210[ 



(2,-2) [148;170] 


(2,-2) [180;190] 


(2,-2) [180; 189] 


(2,-2) [185;192] 
(1,0) [107;110] 
(2,-1) [107;110] 


(2,-2) [175;188] 


- fit within 2 cr; in italics are printed modes without a common range of yot with other frequency peaks 

Model#l - logTeff=4.171; logL/L q=2.368; M = 3.62M©; aov = 0.0 
Model#2 - logrefi=4.157; logL/L0=2.368; M = 3.55M0; aov = 0.0 J7p; 

Model#3 - logTefi=4.171; logL/L0=2.52O; M = 3.84M0; Qov = 0.0 J7p; 

Model#4 - logTefi=4.157; logL/L0=2.52O; M = 3.78M0; Qov = 0.0 J7p; 

Model#5 - logTeff=4.164; logL/L0=2.444; M = 3.64M0; aov = 0.2Hp; 


tailed abundance analysis from high-resolution spectroscopy 
is necessary. 

For some stars with spectroscopic data, viz. HD 25558, 
HD 74195 and HD 182255, we were unable to find solutions 
which reproduce the observed amplitude of the radial veloc¬ 
ity variations to within a la error. Therefore, we accepted 
solutions to within 2a or even 6 a (see Table[^. This can 
be caused by model uncertainties (e. g. model atmospheres, 
chemical composition, opacities, etc. as well as the approach 
of computing stellar pulsations) or observational data (e. g. 
underestimated errors of radial velocity amplitudes or fun¬ 
damental parameters constrained with too low precision). 

Similarly, not always we were able to find a common 
range of the rotational velocity for all modes of a given 
star. The stars in question are: HD 123515, HD 179588, 
HD 208057, and for these stars we followed the case of 
Model^S of our working example HD 21071, i.e., for a given 
frequency we choose a mode with the range of Vrot closest 
to the range of Kot obtained for other frequencies observed 
in the star (they are printed in italic in Tablej^. 

Finally, in some cases the identifications of the mode de¬ 
gree in the zero-rotation approximation and the traditional 
approximation did not agree with each other. An example 
is HD 25558. For this star, both methods indicate that ui 
is a dipole mode, but 1/3 is a mode with £=1, 2 or 3 when 
the zero-rotation approximation is used or £=4, 6 when the 
traditional approximation is used. This is not surprising be¬ 
cause rotation strongly modifies properties of pulsational 
modes of the SPB stars. In particular, rotation changes the 
frequency domain of the unstable modes. This phenomenon 
is especially pronounced for high values of Kot and \m\. 


In the last column of Table|^ we listed for five stars 
spectroscopic determinations of {i , m) taken from the litera¬ 
ture (De Cat et al. 2005). In the case of HD 24587, HD 92287 
and HD 181558, our mode identification is in agreement with 
the angular numbers found by De Cat et al. (2005) but for 
HD 140873 and HD 177863 the results are in contradiction. 
However, if we look at solutions within 3a in for vi in 

HD 177863, the modes (1,0) and (1,4-1) are also possible and 
the last one agrees with the determinations of De Cat et al 
(2005). In the case of vi in HD 140873, we are dealing with a 
similar situation but the agreement is achieved within 4a in 
only additional solution is the mode ( 1 , 4 - 1 ). 
On the other hand, it should be noted that both stars are 
visual and spectroscopic binaries and photometric as well as 
spectroscopic observations can be polluted by the light from 
the companions, and different methods can be sensitive to 
this effect in different ways. Moreover, if ui in HD 177863 is 


tidally excited, then according to Willems & Aerts (2002 1 it 


should be a quadrupole sectoral mode. 


Unfortunately, for HD 21071, which is the best stud¬ 
ied star by us in the sample, there are no determinations 
of angular numbers in the literature from the spectroscopic 
methods and we cannot compare our results. There is only 
the identification of the harmonic degree from the photomet¬ 


ric method. De Cat et al. (20071 found that izi is a dipole 


mode what is in agreement with our result. For U 2 , us and 
1 ^ 4 , they found that these frequencies can be associated with 
dipole or quadrupole modes and these results are also com¬ 
patible with our identifications. The only discrepancy is that 


according to De Cat et al. (20071 U 2 and 1^4 can be also the 
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Table 5. The results of the mode identification for the frequency peaks in all programme stars. The values of the rotational velocity 
and frequency are in km s“^ and d~^, respectively. In the last column (in the section ‘stars with radial velocity variations’), the 
spectroscopically determined values of {i , m) are given if available Cat et al.|2005l l. 


HD 

^ (d-i) 

e{Vrot=0) 

{i,m) 

:1 /..min. ,,max\ 

J i^rot ’ ^rot / 

(e, m) 

Stars with radial velocity variations 

24587 

1.1569 

1 

(1,1) [115;127] {0.6952;0.7677} 


(1,1) 

25558 

0.65265 

1 

(1,0) [174;198] {0.8396;0.9554}2‘^ 




1.93235 

7 

(6,-1) [174;198] {0.8396;0.9554} 




1.17913 

1, 2, 3 

(4,-1) [174;198] {0.8396;0.9554} 

(6,-1) [174;198] {0.8396;0.9554} 


26326 

0.5338 

1, 2 

(1,0) [71;92] {0.3653;0.4733} 




0.7629 

1 

(1,1) [71;92] {0.3653;0.4733} 



53921 

0.6054 

1! 

(2,-1) [17;77] {0.1209;0.5475}! 



74195 

0.35475 

1, 2 

(1,0) [25;43] {0.0986;0.1696}®'^ 

(1,-1) [25;29] {0.0986;0.1144}6'^ 



0.35033 

1, 2 

(1,0) [25;43] {0.0986;0.1696}2‘^ 




0.34630 

1, 2 

(2,-2) [25;43] {0.0986;0.1696} 




0.39864 

1, 2 

(2,-1) [40;43] {0.1578;0.1696} 

(2,-2) [25;40] {0.0986;0.1578} 


74560 

0.6447137 

1 

(1,0) [79;135] {0.4621;0.7897} 

(1,-1) [60;73] {0.3510;0.4270} 



0.3957710 

1 

(1,-1) [60;73] {0.3510;0.4270} 

(1,-1) [79;135] {0.4621;0.7897} 





(2,-1) [115;135] {0.6727;0.7897} 
(2,-2) [110;133] {0.6435;0.7780} 

(2,-2) [60;73] {0.3510;0.4270} 



0.4476463 

1 

(1,-1) [60;73] {0.3510;0.4270} 

(1,-1) [79;135] {0.4621;0.7897} 





(2,-1) [60;73] {0.3510;0.4270} 

(2,-1) [79;135] {0.4621;0.7897} 





(2,-2) [60;73] {0.3510;0.4270} 

(2,-2) [79;90] {0.4621;0.5265} 





(2,-2) [120;135] {0.7020;0.7897} 

(4,-3) [90;91] {0.5265;0.5323} 



0.6359820 

1, 2 

7 




0.8228450 

1 

7 




1.6105673 

1, 2 

7 



85953 

0.2663 

1 

(1,-1) [169;171] {0.5085;0.5145} 




0.2189 

1 

(1,-1) [169;171] {0.5085;0.5145} 

(6,-2) [169;171] {0.5085;0.5145} 


92287 

0.21480 

1 

(1,-1) [122-134] {0.5002;0.5494} 


(1,-1) 

121190 

2.6831 

7 

(1,1) [198-233] {1.7238;2.0285} 




2.6199 

7 

7 




2.4713 

7 

(1,1) [198-233] {1.7238;2.0285} 



123515 

0.68528 

1 

(1,0) [117; 149] {0.5757;0.7332} 




0.65929 

1 

(1,0) [166;194] {0.8169;0.9546} 




0.72585 

1 

(1,0) [166;194] {0.8169;0.9546} 




0.55198 

1, 4 

(1,0) [24;149] {0.1181;0.7332} 



138764 

0.7944 

1, 4 

(1,0) [141;163] {0.9365;!.0826} 

(1,1) [53;64] {0.3520;0.4251} 


140873 

1.1515 

1 

(3,0) [101;102] {0.5921;0.5980} 

(4,0) [80;81] {0.4690;0.4749} 

(1,1) 

177863 

0.84059 

1 

(2,0) [62;63] {0.4074;0.4140} 

(2,0) [68;69] {0.4469;0.4534} 

(1,1) 


0.10108 

7 

(2,-2) [62;63] {0.4074;0.4140} 

(4,-3) [62;63] {0.4074;0.4140} 

(2,-2) [68;69] {0.4469;0.4534} 


181558 

0.80780 

1 

(1,0) [228;242] {1.2699;!.3479} 

(1,1) [77;80] {0.4289;0.4456} 

(1,1) 

182255 

0.97185 

1 

(1,1) [54;75] {0.3921;0.5446}®‘^ 




0.79225 

1 

(1,0) [54;75] {0.3921;0.5446} 




0.62526 

1, 2 

(1,0) [54;65] {0.3921;0.4720} 

(2,-1) [54;75] {0.3921;0.5446} 

(1,-1) [54;75] {0.3921;0.5446} 



1.12780 

1 

(1,0) [54;75] {0.3921;0.5446} 

(2,0) [54;75] {0.3921;0.5446} 

(1,1) [54;75] {0.3921;0.5446} 



1.02884 

1, 2 

(1,0) [54;75] {0.3921;0.5446} 

(1,1) [54;75] {0.3921;0.5446} 




4 

(2,0) [54;75] {0.3921;0.5446} 



208057 

0.89050 

1 

(2,0) [120-122] {0.5935;0.6034} 




0.80213 

1, 2 

(1,0) [138;177] {0.6825;0.8754} 




2.47585 

7 

(4,0) [138;177] {0.6825;0.8754} 



215573 

0.5654 

1, 2 

(1,0) [75;93] {0.4535;0.5623} 

(1,1) [28;33] {0.1693;0.1995} 
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Table 5 — continued 


HD 


e(Vrot = 

=0) {£, m) [l/“-; v;”-] {} 

Stars without radial velocitiy variations 

1976 

0.39946 

1 

i ^5 [143;223] {0.4770;0.7438} 



0.93895 

1 

(1,1) [163;223] {0.5437;0.7438} 

(2,1) [144;163] {0.4803;0.5437} 




(4,0) [143;152] {0.4770;0.5070} 

(6,0) [143;146] {0.4770;0.4870} 


1.20346 

1, 2, 3 

i 7^1 [143;223] {0.4770;0.7438} 


3379 

1.82023 

1, 2 

e^i 



1.59418 

1 



28114 

0.48790 

1 

(1,0) [21;111] {0.1123;0.5935} 

(1,-1) [21;45] {0.1123;0.2406} 




(1,1) [21;41] {0.1123;0.2192} 

(2,0) [21;25] {0.1123;0.1337} 


0.49666 

1 

(1,0) [21;111] {0.1123;0.5935} 

(1,-1) [21;92] {0.1123;0.4919} 




(1,1) [21;41] {0.1123;0.2192} 

(2,0) [21;25] {0.1123;0.1337} 




(2,-1) [21;60] {0.1123;0.3208} 

(3,-1) [73;101] {0.3903;0.5400} 

28475 

0.68369 

1, 2 

(1,0) [30;165] {0.1810;0.9956} 

(1,-1) [30;74] {0.1810;0.4465} 




(1,1) [30;55] {0.1810;0.3319} 

(2,-1) [30;63] {0.1810;0.3802} 




(2,0) [30;36] {0.1810;0.2172} 

(2,0) [38;50] {0.2293;0.3017} 




(3,-1) [41;42] {0.2474;0.2534} 

(4,-1) [85;206] {0.5129;1.2430} 


0.40893 

1 

(1,-1) [30;203] {0.1810;1.2249} 

(2,-1) [105;156] {0.6336;0.9413} 




(2,-2) [30;56] {0.1810;0.3379} 

(2,-2) [110;116] {0.6638;0.7000} 




(4,-2) [50;54] {0.3017;0.3258} 

(r,-l) [165;206] {0.9956;!.2430} 

(4,-3) [31;36] {0.1871;0.2172} 

34798 

0.78350 

1 

(1,0) [31;157] {0.1878;0.9513} 

(1,-1) [30;49] {0.1818;0.2969} 




(1,1) [30;65] {0.1818;0.3939} 

(2,0) [40;70] {0.2424;0.4242} 

(2,0) [30;34] {0.1818;0.2060} 


0.85231 

1, 2, 4 

? 



0.97432 

1 

(1,0) [31;157] {0.1878;0.9513} 

(1,1) [30;100] {0.1818;0.6059} 




(2,0) [30;35] {0.1818;0.2121} 

(2,0) [39;120] {0.2363;0.7271} 




(2,1) [30;50] {0.1818;0.3030} 

(2,2) [30;30] {0.1818;0.1818} 




(3,0) [44;70] {0.2666;0.4242} 

(4,-1) [89;157] {0.5393;0.9513} 




(5,-1) [129;157] {0.7817;0.9513} 

(6,-6) [30;30] {0.1818;0.1818} 


0.67795 

1, 2 

(1,0) [30;145] {0.1818;0.8786} 

(1,-1) [30;101] {0.1818;0.6120} 




(1,1) [30;45] {0.1818;0.2727} 

(2,0) [30;35] {0.1818;0.2121} 




(2,-1) [30;75] {0.1818;0.4545} 
(4,-2)[30;30] {0.1818;0.1818} 

(3,-1) [36;157] {0.2181;0.9513} 

37151 

1.2431665 

1! 

(1,0) [20;100] {0.2195;1.0974} 

(1,1) [20;45] {0.2195;0.4938} 

(1,-1) [20;140] {0.2195;1.5363} 


1.1801861 

1 

(1,0) [20;85] {0.2195;0.9328} 

(1,1) [20;40] {0.2195;0.4390} 

(1,-1) [20;140] {0.2195;1.5363} 


1.0426632 

1 

(1,0) [20;60] {0.2195;0.6584} 

(1,1) [20;25] {0.2195;0.2743} 

(1,-1) [20;140] {0.2195;1.5363} 


1.1052525 

1 

(1,0) [20;70] {0.2195;0.7682} 

(1,-1) [20;140] {0.2195;1.5363} 




(1,1) [20;30] {0.2195;0.3292} 

(2,-2) [30;35] {0.3292;0.3841} 


1.4520631 

1, 2 

(1,0) [20;140] {0.2195;1.5363} 

(1,-1) [20;70] {0.2195;0.7682} 




(1,1) [20;65] {0.2195;0.7133} 

(2,-2) [20;68] {0.2195;0.7462} 

(2,-1) [20;105] {0.2195;1.1522} 

45284 

1.23852 

1 

(1,0) [53;272] {0.4083;2.0954} 

(1,1) [52;102] {0.4006;0.7858} 




(2,0) [52;60] {0.4006;0.4622} 

(2,0) [67;117] {0.5161;0.9013} 




(3,0) [62;67] {0.4776;0.5161} 

(5,-1) [119;152] {0.9167;!.1709} 

(4,-1) [103;272] {0.7935;2.0954} 


1.12753 

1 

(1,0) [53;272] {0.4083;2.0954} 

(1,-1) [52;52] {0.4006;0.4006} 




(1,1) [52;87] {0.4006;0.6702} 

(2,0) [52;60] {0.4006;0.4622} 




(2,0) [68;97] {0.5238;0.7472} 

(4 ,-l) [106;192] {0.8166;1.4791} 


1.50605 

? 

7 
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Table 5 — continued 


HD V (d-i) 

^(Rrot=0) 

(£,m) 

1 / ,,min.-.max \ 

1 l^rot 5 ^rot / 

143309 0.59830 

1 

(1,0) [10;121] {0.0611;0.7392} 

(1,1) [10;45] {0.0611;0.2749} 

(2,0) [10;11] {0.0611;0.0672} 

(2,-2) [10;12] {0.0611;0.0733} 

(1,-1) [10;19] {0.0611;0.1161} 
(2,-1) [10;30] {0.0611;0.1833} 
(2,0) [14;30] {0.0855;0.1833} 

0.60213 

1, 2 

(1,0) [10;121] {0.0611;0.7392} 

(1,1) [10;45] {0.0611;0.2749} 

(2,0) [13;30] {0.0794;0.1833} 

(2,-2) [10;14] {0.0611;0.0855} 

(1,-1) [10;22] {0.0611;0.1344} 
(2,0) [10;11] {0.0611;0.0672} 
(2,-1) [10;40] {0.0611;0.2444} 

0.71414 

1, 2 

(1,0) [10;121] {0.0611;0.7392} 

(1.1) [10;65] {0.0611;0.3971} 

(2,0) [13;60] {0.0794;0.3666} 

(2.1) [10;25] {0.0611;0.1527} 

(2.2) [10;15] {0.0611;0.0916} 

(1,-1) [10;23] {0.0611;0.1405} 
(2,0) [10;11] {0.0611;0.0672} 
(2,-1) [10;42] {0.0611;0.2566} 
(2,-2) [10;15] {0.0611;0.0916} 

1.20066 

1 

(2,0) [74;121] {0.4521;0.7392} 

(3,0) [63;115] {0.3849;0.7026} 

(4,1) [21;35] {0.1283;0.2138} 

(4,4) [10;10] {0.0611;0.0611} 

(2.1) [29;90] {0.1772;0.5498} 
(4,0)[71;80] {0.4338;0.4887} 

(4.2) [14;20] {0.0855;0.1222} 

0.59133 

1, 2 

(1,0) [10;120] {0.0611;0.7331} 

(1,1) [10;45] {0.0611;0.2749} 

(2,0) [13;30] {0.0794;0.1833} 

(2,-2) [10;22] {0.0611;0.1344} 

(1,-1) [10;54] {0.0611;0.3299} 
(2,0) [10;12] {0.0611;0.0733} 
(2,-1) [10;121] {0.0611;0.7392} 
(3,-1) [35;95] {0.2138;0.5804} 


160124 0.52014 1 


0.52096 1 


0.52055 1 


0.52137 ? 

0.52197 1 


0.70259 1 


1.0422 1, 2, 3, 4 

0.31366 1, 2 


(1,0) [11;78] {0.0439;0.3113} 

(1.1) [74;78] {0.2954;0.3113} 
(2,0) [13;78] {0.0519;0.3113} 

(2.2) [8;14] {0.0319;0.0559} 
(1,0) [10;78] {0.0399;0.3113} 

(1.1) [63;78] {0.2514;0.3113} 
(2,0) [12;78] {0.0479;0.3113} 

(2.2) [8;17] {0.0319;0.0679} 
(1,0) [12;78] {0.0479;0.3113} 

(1.1) [76;78] {0.3033;0.3113} 

(2.1) [8;38] {0.0319;0.1517} 
(3,0) [42;58] {0.1676;0.2315} 
? 

(1,0) [11;78] {0.0439;0.3113} 

(1.1) [74;78] {0.2954;0.3113} 
(2,0) [13;78] {0.0519;0.3113} 

(2.2) [8;14] {0.0319;0.0559} 
(2,0) [8;8] {0.0319;0.0319} 
(2,-1) [9;26] {0.0359;0.1038} 
(2,-2) [8;10] {0.0319;0.0399} 

(1,0) [8;78] {0.0319;0.3113} 
(1,1) [8;23] {0.0319;0.0918} 
(2,-2) [8;19] {0.0319;0.0758} 


(1.1) [8;11] {0.0319;0.0439} 
(2,0) [8;8] {0.0319;0.0319} 

(2.1) [8;38] {0.0319;0.1517} 
(3,0) [40;58] {0.1596;0.2315} 

(1.1) [8;14] {0.0319;0.0559} 
(2,0) [8;8] {0.0319;0.0319} 

(2.1) [8;38] {0.0319;0.1517} 
(3,0) [36;58] {0.1437;0.2315} 

(1.1) [8;10] {0.0319;0.0399} 
(2,0) [13;78] {0.0519;0.3113} 

(2.2) [8;14] {0.0319;0.0559} 


(1.1) [8;11] {0.0319;0.0439} 
(2,0) [8;8] {0.0319;0.0319} 

(2.1) [8;38] {0.0319;0.1517} 
(3,0) [42;58] {0.1676;0.2315} 
(2,0) [13;30] {0.0519;0.1197} 
(2,1) [59;78] {0.2355;0.3113} 
(4,0) [64;78] {0.2554;0.3113} 

(1,-1) [8;45] {0.0319;0.1796} 
(2,-1) [8;67] {0.0319;0.2674} 


160762 

0.28671 

7 

(1,-1) [98;105] {0.4414;0.4729} 


179588 

0.85654 

1, 4 

(1,0) [125;205] {0.6481;1.0628} 



2.04263 

1!, 2! 

(2,1) [177; 190] {0.9177;0.9851}!! 



2.19989 

1 

(2,1) [203;205J {1.0525; 1.0628} 



1.83359 

1 

(2,1) [135; 165] {0.6999;0.8555} 


191295 

0.71908 

1 

(1,0) [17;128] {0.1033;0.7774} 

(1,-1) [16;32] {0.0972;0.1944} 




(1.1) [16;71] {0.0972;0.4312} 

(2.1) [16;26] {0.0972;0.1579} 

(2,0) [37;56] {0.2247;0.3401} 


0.49011 

1, 2 

(1,0) [16;76] {0.0972;0.4616} 

(1,-1) [16;128] {0.0972;0.7774} 




(1,1) [16;31] {0.0972;0.1883} 

(2,-2) [16;35] {0.0972;0.2126} 

(2,-1) [16;81] {0.0972;0.4920} 


0.29302 

1 

(1,-1) [21;128] {0.1275;0.7774} 

(2,-2) [31;49] {0.1883;0.2976} 
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Table 5 — continued 


HD V (d-i) 

^(Uot=0) 

{£,m) [KS:“;v;or] 

206540 0.72002 

1 

(1,0) [17;149] {0.0932;0.8165} 
(1,-1) [15;21] {0.0822;0.1151} 
(2,0) [59;85] {0.3233;0.4658} 

(2,2) [15;25] {0.0822;0.1370} 

(1,0) [160;177] {0.8768;0.9699} 

(1.1) [15;85] {0.0822;0.4658} 

(2.1) [21;40] {0.1151;0.2192} 

0.62125 

1, 2, 4 

(1,0) [15;149] {0.0822;;0.8165} 
(1,-1) [15;35] {0.0822;0.1918} 
(2,0) [15;17] {0.0822;0.0932} 
(2,-1) [15;76] {0.0822;0.4165} 
(2,-2) [15;21] {0.0822;0.1151} 

(1,0) [160;177] {0.8768;0.9699} 

(1.1) [15;65] {0.0822;0.3562} 

(2,0) [20;60] {0.1096;0.3288} 

(2.1) [15;25] {0.0822;0.1370} 

(2.2) [15;15] {0.0822;0.0822} 

0.38271 

1, 2 

(1,0) [15;55] {0.0822;0.3014} 

(1,1) [15;15] {0.0822;0.0822} 
(2,-2) [15;36] {0.0822;0.1973} 

(1,-1) [15;149] {0.0822;0.8165} 
(2,-1) [25;75] {0.1370;0.4110} 
(r,-l) [160;177] {0.8768;0.9699} 


- radial velocity within 2(7; - radial velocity within 6(7; 

! — accepted ^ 3.3; !! - accepted 'p Si 4.6; 

? - no constraints; in italics are printed modes without a common range of l4ot with other frequency peaks. 


£ = 4 mode. In fact, this solution was found also by us but 
we rejected it because this mode was stable. 


4 SUMMARY 


We obtained a mode identification or at least constraints 
on angular numbers (t, m) for 31 SPB stars. Simultane¬ 
ously, we derived bounds on the rotational velocities of the 
stars we examined. It should be noted, that the results de¬ 
livered by the presented method depend on the model pa¬ 
rameters. It would be very interesting to verify them on the 
SPB stars for which the inclination angle and/or rotation 
period are known from independent determinations. We see 
such prospects in the forthcoming two-colour photometric 
data from the BRITE satellites and follow-up spectroscopic 
observations. 

To summarize our results, we constructed histograms 
showing the distributions of the mode degree, £, (Fig.[^, 
azimuthal order, m (Fig.j^, the intrinsic mode amplitude, 
£ (Fig.[^, the values of Kot and Kotsini (Fig.[^, and the 
inclination angle, i (Fig.[^. Due to the fact that for some 
frequencies the identification was not unambiguous, we shall 
now explain how we have counted the modes with a given 
value of the mode degree or azimuthal order. In the case of 
the modes for which we found j combinations of (£, m) and 
£ ox m occurred i times, we counted such angular numbers 
i/j times. The histograms for Kot and i were constructed in 
a similar way. 

When we neglect all effects of rotation, most of the 
modes appeared to be dipole modes (see Fig.|^. Some modes 
were quadrupole modes and very few had I equal to 3 or 4. 
None of the frequencies had f = 5 or 6. If the effects of 
rotation were included, the number of the dipole modes de¬ 
creased significantly and number of the quadrupole modes 
was approximately doubled. We also obtained slightly larger 
number of the modes with I higher than 2. Despite this, only 
very few modes had One can see also that there were 

fewer modes with £ = 3 and 5 than modes with ^ = 4 and 6. 
This is because the odd degree modes have lower visibility 
(see Daszyhska-Daszkiewicz et al.|2002 20151. 



1 2 3 4 5 6 

£ 


Figure 5. The number of the modes with different values of the 
harmonic degree, I, obtained with (black filled bars) and without 
(hatched bars) the effects of rotation taken into account. 


From Fig.j^one can see that zonal modes are quite nu¬ 
merous amongst our sample of the 31 stars. We have to 
admit that we do not have any explanation of that fact. 
Maybe this is the case or it can result from model uncer¬ 
tainties and/or limitations of the method we applied. The 
number of sectoral retrograde and prograde modes is about 
the same and if we count them all, the number of the sec¬ 
toral modes becomes only slightly smaller than the number 
of the zonal modes (cf. the insert of Fig.|^. Some modes 
with higher values of m also happen. 

In the method we used, the values of the intrinsic mode 
amplitude, e, are derived from observations (equation 3) and 
used to constrain mode identification. The histogram for e is 
shown in Fig.[^ An interesting fact is that most values of e 
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Figure 6. The number of the modes with different values of the 
azimuthal order, m. The insert shows the histogram of the number 
of the zonal (z), sectoral retrograde (sr), sectoral prograde (sp) 
and tesseral (t) modes. 


Figure 8. Distribution of the values of Frotsinf (hatched bars) 
and the values of Vl-ot (black solid bars) we obtained simultane¬ 
ously with the identification of {I, m) for the studied stars. 



0 0.002 0.004 0.006 0.008 


|el 

Figure 7. Distribution of the values of e for the frequency peaks 
of the studied stars. 


10 



0 15 30 45 60 75 90 


Figure 9. Distribution of the inclination angle, i, for the stars of 
Table 1. 


are much smaller than the adopted upper limit, £max = 0.01. 
A similar result was obtained for the frequency peaks of 
/t Eri by Daszynska-Daszkiewicz et al. (20151. In the case 
of our sample of stars, above 60 per cent of the frequency 
peaks have e below 0.005. If we limit the sample to the stars 
with the radial velocity measurements, the number of the 
frequency peaks with e < 0.005 increases to above 80 per 
cent. 

In Fig.|^ we present a histogram showing the distri¬ 


bution of the values of the rotational velocity we obtained 
and the projected rotational velocity. One can see that a 
substantial number of the programme stars have Kot sin i 
below 40 km s“^. 

One can observe a quick decrease of the number of stars 
with large values of Kotsini, and there are no stars with 
Vrot sini higher than 160 km s“^. The view is quite different 
when we consider the ‘true’ value of the rotational velocity. 
In the first place, stars rotating below 40 km s“^ are not 
most numerous and there are stars with rotational velocities 
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from the whole range of Kot • The number of stars with K-ot 
in the range 40-80 km s“^ is greater by a factor of 2 than 
that of stars with the smallest rotational velocities. Then 
the number of stars is slowly decreasing with the increasing 
value of Kot- 

The distribution of the projected rotation velocities and 
rotation velocities of B-type stars can be found in, e. g. 


Huang & Gies (20061. Although in the cited paper it is 
shown that the distribution of Kot sin i for the field and 
cluster stars differ slightly, it is striking that the values of 
Vrot sin i in our sample are lower and the number of stars 
quickly decreases with increasing Kot sin i. [Huang fc Gies| 
(20061 also showed that most numerous among B-type stars 
are those rotating with the equatorial velocity equal to about 
200 km s“^. In our sample, the most numerous are stars ro¬ 
tating in the range 40-80 km s“'^. Such shift in distribution 
can suggest that SPB stars are indeed slower rotators than 
the other B-type stars (e.g. Balona 20101. As opposed to 


Huang & Gies (20061, Levato & Grosso (20131 found that 


the distribution of K-ot sin i as well as Kot for B-type stars is 
bimodal what is not seen in our sample. However, we want to 
emphasize that our sample is too small to make any certain 
conclusions. 

Finally, in Fig.is presented a distribution of the incli¬ 
nation angle, i, we derived for our sample of stars. One can 
see that about half of the stars are observed at rather low 
inclination angle i < 30°. As shown by Abt (20011, the stel¬ 
lar rotational axes should be distributed randomly so one 
would expect rather the flat distribution of i but again we 
have to remind that we are dealing with a small sample and 
the shape of histogram can be accidental. The other reason 
could be that we identified so many zonal modes which are 
best seen from the pole-on direction (cf. Fig.[^. The other 
explanation is that the assumption of uniform rotation is 
incorrect and the higher values of Vj-ot (implying the low 
values of i) correspond to the deeper layers which are better 
probed by high-order g modes. Some support for this hy¬ 
pothesis is the fact that there is a comparable population of 
the sectoral modes which are least cancelled when we look at 
the equator and there is no increase in the number of stars 
with the inclination angle close to 90°. This is an interesting 
explanation but it demands more detailed studies which go 
beyond the scope of the present paper. 

The unambiguity of mode identification strongly de¬ 
pends on the quality of the photometric and/or spectro¬ 
scopic data. The quality of the data is especially important 
when the effects of rotation are included. This is due to the 
fact that contrary to the case of the zero-rotation approxi¬ 
mation when we determined only one parameter, I, in the 
case with rotation included we have to obtain additionally 
the azimuthal order, m, and rotational velocity, Kot. The 
information contained in the amplitude of the radial veloc¬ 
ity variations helps to reduce significantly the number of 
possible solutions. Thus, combined photometric and spec¬ 
troscopic data are highly desirable. 

In the case of HD 123515, HD 179588 and HD 208057 
as well as in Model/!^3 and Model/fi5 of HD21071, we failed 
Ending a common range of the rotational velocity for all ob¬ 
served frequencies. This can be explained in several ways. 
First, the accuracy of our pulsational calculations may be 
insufficient and the theory needs to be improved. Secondly, 
the formal errors of the amplitudes and phases of the light 


and radial velocity variations were underestimated. Thirdly, 
the astrophysical parameters we used, logTefi and logL/LQ, 
were not possible to be constrained with a high precision. 
It is also possible that the adopted initial abundance of hy¬ 
drogen, metallicity or overshooting from the convective core 
made the starting models inappropriate. And again, it is 
possible that our assumption of rigid rotation is not sat¬ 
isfied, so that different K-ot ranges for different frequency 
peaks observed in a star are manifestations of differential 
rotation. Perhaps all above-mentioned factors play a role 
simultaneously. 

Finally, close companions can disrupt photometric and 
spectroscopic measurements and thus the mode identifica¬ 
tion. Two of three stars with the high solutions, HD 53921 
and HD 179588, are visual binary or visual and spectroscopic 
binary, respectively. A high value of can also suggest non- 
pulsational originate of variability. 

Nevertheless, we see prospects for seismic modelling of 
the SPB stars for which the best constraints on the angu¬ 
lar numbers, {i, m) and the rotational velocity, V(ot, were 
obtained. Only in such cases one can try to attribute the 
radial order, n, to the observed frequency peaks. In par¬ 
ticular, HD 21071 is a promising candidate. Our results of 
seismic modelling of this star will be presented in a separate 
paper (Szewczuk & Daszyriska-Daszkiewicz, in preparation). 
Preliminary results has been already presented injSzewczukj 


& Daszyhska-Daszkiewicz (20151. 
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